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I .  INTRODUCTION 

The  parametric  receiver  utilizes  the  nonlinearity  of  acoustic  wave 
propagation  in  water  to  synthesize  a  virtual  array  between  two  trans¬ 
ducers,  the  pump  and  hydrophone  (see  Appendix  A  for  details).  This  vir¬ 
tual  array  has  the  directivity  characteristics  of  a  conventional  end-fired 
array  of  the  same  length  as  the  pump-hydrophone  separation,  but  with  a 
considerable  reduction  in  the  number  of  transducers  required.  The  main 
advantage  of  using  a  parametric  receiver  in  a  submarine  application, 
therefore,  is  that  a  narrow,  conical  receiving  beam  can  be  obtained  with 
minimum  hardware  in  the  water.  This  receiving  beam  has  a  front-to-back 
ratio  and  a  vertical  directivity  that  discriminate  against  noise  at  the 
mid-frequencies.  Also,  since  the  virtual  array  is  synthesized  in  the 
water,  it  is  reformed  (or  stabilized)  at  the  speed  of  sound  during  maneu¬ 
vers  and  heading  changes.  Because  there  are  no  transducer  elements  between 
pump  and  hydrophone,  the  parametric  receiver  may  be  less  sensitive  to  flow 
noise  as  boat  speed  increases.  These  characteristics  of  the  parametric 
receiver  contribute  to  its  attractiveness  in  submarine  passive  sonar 
applications.  However,  there  are  possible  performance  limitations  of 
mobile  parametric  receivers  that  need  to  be  considered. 

Because  the  parametric  receiver  effectively  forms  an  array  in  the 
medium  between  pump  and  hydrophone,  its  operation  is  sensitive  to  the 
state  of  the  medium  in  this  region.  The  presence  of  air  bubbles,  pieces 
of  hardware,  or  significant  temperature  or  flow  velocity  gradients  will 
modify  the  synthesized  array  and  alter  its  performance.  In  particular, 
at  moderate  or  high  boat  speeds,  the  turbulence  produced  in  the  wake  of 
the  pump  transducer  may  have  several  effects  on  the  performance  of  the 
parametric  receiver.  The  turbulence  will  act  as  a  volume  distributed, 
low  frequency  acoustic  source  that  will  increase  the  noise  level  of  the 
parametric  receiver.  Turbulence  will  also  scatter  the  acoustic  waves 
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propagating  between  pump  and  hydrophone,  thereby  producing  random 
amplitude  and  phase  fluctuations  in  the  detected  acoustic  signals.  These 
effects  of  turbulence  need  to  be  Investigated  in  order  to  determine  the 
viability  of  the  parametric  receiver  as  a  submarine  passive  sonar. 

The  objective  of  the  present  study  is  to  analytically  determine  the 
effects  of  turbulence  on  the  performance  of  the  mobile  parametric  receiver, 
with  emphasis  on  the  effects  of  turbulent  scattering.  The  study  has  been 
organized  into  three  general  areas.  (1)  A  literature  survey  has  been  con¬ 
ducted  in  the  areas  of  parametric  reception,  turbulence,  and  wave  propa¬ 
gation  in  inhomogeneous  media.  Summaries  of  material  from  these  areas 
that  pertain  to  the  study  are  presented  in  Appendices  A,  B,  and  C. 

(2)  Theoretical  expressions  have  been  developed  for  the  acoustic  amplitude 
and  phase  fluctuations  caused  by  turbulent  scattering.  A  summary  and 
discussion  of  this  theory  is  presented  in  Section  II  of  this  report,  with 
details  of  the  analysis  given  in  Appendices  D  and  E.  (3)  Approximate 
numerical  estimates  of  the  effects  of  turbulence  on  parametric  reception 
have  been  made  using  the  developed  theory.  These  estimates  are  presented 
and  discussed  in  Section  III  of  the  report.  Finally,  in  Section  IV,  con¬ 
clusions  are  drawn  regarding  the  limitations  imposed  by  turbulence  on 
mobile  parametric  reception,  and  future  work  is  discussed. 
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II.  THEORETICAL  RESULTS 


A.  Introduction 


Suppose  a  parametric  receiver  consisting  of  pump  transducer  and 
hydrophone  is  placed  on  a  vessel  moving  from  right  to  left  through  the 
medium.  If  the  transducers  are  placed  clear  of  the  boundary  layer  tur¬ 
bulence  produced  by  the  vessel's  hull,  then  the  situation  can  be  modeled 
by  considering  two  stationary  transducers  in  a  fluid  that  flows  from  left 
to  right,  as  shown  in  Fig.  1.  The  flow  will  separate  in  the  region  of 
the  pump  transducer,  and  vortices  will  form.  At  some  distance  downstream 
from  the  transducer,  the  vortices  will  decompose  into  the  random  velocity 
field  that  characterizes  turbulence.  The  dimensions  and  intensity  of  the 
turbulence  are  dependent  upon  the  flow  velocity  and  upon  the  geometry  of 
the  rigid  boundaries  associated  with  the  flow. 

As  shown  in  Appendix  B,  the  eddies  associated  with  turbulent  flow 
may  be  treated  as  "patches"  of  variable  refractive  index.  These  patches, 
or  inhomogeneities,  will  scatter  the  acoustic  waves  propagating  in  the 
interaction  region  between  pump  transducer  and  hydrophone.  As  a  result 
of  scattering,  amplitude  and  phase  fluctuations  in  both  pump  and  sideband 
waves  will  occur  (see  Appendix  C) . 

The  fluctuations  in  the  sideband  waves  are  a  source  of  noise  to  the 
parametric  receiver  that  can  act  to  degrade  its  performance  in  the  detec¬ 
tion  of  acoustic  signals.  The  purpose  of  the  theoretical  work  discussed 
in  this  section  is  to  obtain  expressions  for  the  amplitude  and  phase 
fluctuations  in  the  sideband  waves  in  terms  of  parameters  of  the 
turbulence . 


B.  Summary  of  Theoretical  Analysis 


1.  Analysis  for  Weak  Turbulence 


A  detailed  analysis  of  the  effects  of  weak  turbulence  on  the 
performance  of  a  mobile  parametric  receiver  is  given  in  Appendix  D.  In 
this  analysis,  the  acoustic  waves  are  modeled  as  follows.  The  pump  wave 
is  assumed  to  be  spherically  spreading,  but  confined  to  a  narrow  beam  by 
the  directivity  of  the  pump  transducer.  There  are  fluctuations  and 
Sp  in  the  amplitude  and  phase,  respectively,  of  the  pump  wave.  For  the 
case  of  weak  turbulence,  these  fluctuations  are  assumed  to  be  small  (i.e., 
less  than  10%  of  the  mean  amplitude  or  phase).  The  signal  wave  is  assumed 
to  be  planar  and,  because  of  its  relatively  low  frequency,  is  assumed  to 
have  a  negligible  level  of  fluctuations  due  to  the  turbulence  in  the 
interaction  region. 


The  interaction  of  these  two  first-order  waves  (the  pump  and 
signal  waves)  produces  an  array  of  virtual  sources  in  the  region  between 
pump  and  hydrophone.  The  second-order  pressure  radiated  from  each  vir¬ 
tual  source  will  have  amplitude  and  phase  fluctuations,  B+  and  S+, 
respectively,  due  to  scattering  caused  by  turbulence.  The  second-order 
pressure  p+  detected  by  the  hydrophone  of  the  parametric  receiver  will 
therefore  depend  upon  the  fluctuation  terms  B^,  S^,  B+,  and  S+. 

The  mean  squared  amplitude  and  phase  fluctuations  in  the 
pressure  p+  are  shown  in  Appendix  D  to  be 


<b{L>  =  0.1415  C2  k7/6L11/6 
PR  n  p 


<S2  >  =  0.8889  /n  <uZ>  kZ  aL  for  D>>1,  and 
r  R  P 


2.  ,  2 


<S2r>  =  1.778  Jn  <u2>  k2  aL  for  D«1 


(1) 

(2) 
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where 


C  is  the  turbulence  structure  constant, 
n 

k  is  the  acoustic  wave  number  at  frequency  u>  , 

P  P 

L  is  the  pump-hydrophone  separation, 

2 

<U  >  is  the  mean  squared  refractive  index  variation, 

a  is  the  mean  correlation  distance  of  the  refractive  index,  and 

„  4L 
D  =  — 

ka“ 


Equations  (1)  and  (2)  are  expressions  for  fluctuations  in  the  detected 

pressure  p+  in  terms  of  the  acoustic  wave  number,  the  array  length  L,  and 

2 

the  turbulence  parameters  C  ,  <p  >,  and  a.  These  results  indicate  that 

the  fluctuations  increase  with  the  intensity  of  the  turbulence  in  the 

interaction  region.  This  is  because  the  amplitude  and  phase  fluctuations 
2  2 

increase  with  and  <u  >,  respectively,  and  both  of  these  parameters  are 
related  to  turbulent  intensity.  Also  from  Eqs.  (1)  and  (2),  it  can  be 
seen  that  the  fluctuations  increase  with  the  separation  L  of  the  pump  and 
hydrophone.  This  is  reasonable,  because  as  L  increases  the  number  of 
scatterers  that  lie  in  the  paths  of  the  propagating  waves  is  increased. 

Equations  (1)  and  (2)  were  obtained  using  the  assumption  that 
complete  longitudinal  correlation  of  fluctuations  exists  in  the  inter¬ 
action  region  of  the  parametric  receiver.  This  assumption  has  the  effect 
of  making  the  equations  apply  for  the  "worst  case",  and  thus  give  maximum 
values  of  fluctuations.  In  applications  where  the  pump-hydrophone  separa¬ 
tion  is  greater  than  the  correlation  distance  of  the  pump  wave  fluctuations, 
2 

and  <S  >  will  be  less  than  values  obtained  from  Eqs.  (1)  and  (2). 

"K  rK  9 

General  expressions  for  calculating  <B_>  and  <S“  >  are  given  in 

PK  PR 

Appendix  D  [Eqs.  (D-24)  and  (D-30) ] . 

The  results  are  also  based  on  the  assumption  that  the 
fluctuations  in  the  pump  and  interaction  frequency  waves  remain  small. 

This  amounts  to  assuming  that  the  turbulence  is  sufficientl”  weak  for 
given  array  lengths  and  acoustic  wa\e  numbers  that  small  perturbation 
theory  is  applicable  to  the  problem.  For  stronger  turbulence,  where  the 


methods  of  small  perturbation  theory  no  longer  apply.  It  Is  useful  to 
make  a  simpler  theoretical  model  than  that  used  in  Appendix  D. 


2.  Analysis  for  Strong  Turbulence 


It  is  shown  in  Appendix  E  that  the  fluctuations  in  the  sideband 
pressure  p+  are  approximately  equal  to  corresponding  fluctuations  in  the 
pump  wave.  In  other  words. 


<BPR> 


<b2> 

p 


and 


> 


<S2> 


As  an  approximate  solution  to  the  problem  of  strong  turbulence,  smooth 
perturbation  theory  for  linear  waves  may  be  used  to  obtain  the  following 
results  (see  Appendix  E  for  details) : 


> 


0.13  C2  k7/6  L11/6 
n  p 


(3) 


and 


<S2  >  =  0.50  <y2>  k2  aL  .  (4) 

rK  P 

It  can  be  seen  that  these  results  are  approximately  equal  to  the  expres¬ 
sions  in  Eqs.  (1)  and  (2),  the  difference  being  in  the  numerical  constants 
This  difference  is  more  severe  for  the  phase  fluctuations  than  for  the 
amplitude  fluctuations.  However,  Eqs.  (3)  and  (4)  are  valid  for  rms  fluc¬ 
tuation  levels  up  to  about  50%,  whereas  Eqs.  (1)  and  (2)  are  valid  only 
up  to  rms  fluctuation  levels  of  about  10%.  The  simple  "strong  turbulence" 
model  described  in  Appendix  E  thus  has  the  effect  of  extending  the  range 

of  validity  of  the  theoretical  results  obtained  in  Appendix  D,  although 
2 

for  <Sp„>  the  approximation  is  somewhat  crude. 


III.  THEORETICAL  EXAMPLES 


The  results  summarized  in  Section  II  can  be  used  to  predict  rms 

2  1/2  2  1/2 

levels  of  fluctuations,  <BpR>  and  <SpR>  ,  for  a  mobile  parametric 
receiver.  As  the  theoretical  results  depend  strongly  upon  parameters  of 
turbulence,  some  estimates  for  these  parameters  need  to  be  made.  It  is 
shown  in  Appendix  F  that  the  pertinent  parameters  may  be  calculated  as 
follows.  The  structure  constant  is  given  by 

C2  =  8.487  x  10-8  v2  L-2/3  ,  (5) 

n  o 


where 

v  is  the  velocity  of  the  vessel  in  m/sec,  and 

L  is  the  outer  scale  of  turbulence,  determined  by  dimensions  of  the 
o 

flow  around  the  pump  transducer  housing. 

The  mean  square  refractive  index  variations  are 

<u2>  =  ^0.3162  v/cQj2  ,  (6) 


where  c  is  the  mean  sound  speed  in  m/sec.  The  mean  correlation  distance 
o 

associated  with  the  refractive  index  variations  is  given  by 

a  =  4.737  x  10'2  AT  .  (7) 

o 

21/2  2  1/2 

Predictions  of  <B  >  and  <S„  >  can  be  made  by  substituting  Eqs.  (5)- 
PR  PR 

(7)  into  Eqs.  (1)  and  (4).  Results  of  such  predictions  appear  in  Figs.  2, 
3,  and  4.  Note  that  "strong  turbulence  theory"  is  used  for  the  phase 
fluctuations  because  the  fluctuation  levels  exceed  0.1  rad  for  most  cases 
considered. 
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FIGURE  3 
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In  Fig.  2,  values  of  rms  amplitude  fluctuations  <BpR>  are  shown 

as  a  function  of  array  length  L  and  pump  frequency  f^.  These  values  are 

calculated  assuming  v  ■  5  kt  and  Lq  ■  0.25  m  (here  Lq  is  approximated  as 

the  diameter  of  the  wake  behind  the  pump  transducer) .  It  can  be  seen 

that,  even  at  the  low  speed  of  5  kt,  significant  levels  of  rms  amplitude 

fluctuations  are  predicted  for  array  lengths  greater  than  10  m  and  pump 

frequencies  greater  than  100  kHz.  For  example,  with  L  *  20  m  and 

2  1/2 

f  ■  250  kHz,  the  predicted  value  of  <B  >  is  approximately  0.4.  This 

P  rR 

means  that  the  "output"  of  the  parametric  receiver  (i.e.,  the  sideband 

pressure  detected  by  the  hydrophone)  will  vary  in  amplitude  by  an  rms 

amount  of  40%  of  its  mean  value.  As  shown  in  the  figure,  this  value  of 

fluctuation  increases  with  increasing  pump  frequency. 

2  1/2 

Similar  results  are  shown  in  Fig.  3  for  <B  >  as  a  function  of 

2  R 

array  length  L  and  structure  constant  C^.  Shown  in  parentheses  are 
approximate  values  of  boat  speed  to  which  the  structure  constants  cor¬ 
respond.  The  pump  frequency  is  assumed  to  be  500  kHz  and  Lq  is  again 
0.25  m.  It  can  be  seen  that  increasing  boat  speed  produces  a  significant 
Increase  in  the  level  of  amplitude  fluctuation  for  a  given  array  length. 

For  an  array  length  of  20  m,  the  level  of  amplitude  fluctuation  exceeds 

2  1/2 

the  range  of  validity  of  the  theory  (i.e.,  <B_D>  exceeds  0.5)  for 
2 

most  values  of  C  shown. 

n 

Phase  fluctuations  as  a  function  of  array  length  and  pump  frequency 

2  1/2 

are  shown  in  Fig.  4.  These  values  of  rms  phase  fluctuations  <S  >  , 

2  -1 

in  radians,  are  calculated  for  <p  >  =  2.94  *  10  ',  which  corresponds  to  a 
boat  speed  of  approximately  5  kt,  and  for  Lq  =  0.25  m.  The  fluctuation 
levels  in  the  figure  are,  for  an  array  length  of  20  m,  generally  in  excess 
of  0.1.  Thus  it  can  be  seen  from  Fig.  4  that,  for  pump  frequencies 
greater  than  250  kHz  and  array  lengths  greater  than  20  m,  the  phase  fluc¬ 
tuations  will  have  significant  levels  (in  excess  of  0.35  radians).  The 
levels  shown  in  Fig.  4  may  be  expected  to  increase  with  boat  speed  In  a 

manner  similar  to  that  demonstrated  by  the  amplitude  fluctuations.  For 

2  1/2 

a  pump  frequency  of  500  kHz,  <SpR>  will  exceed  0.5  rad  for  speeds 
greater  than  5  kt,  even  at  short  array  lengths. 
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IV.  CONCLUSIONS  AND  FUTURE  WORK 

The  objective  of  the  investigation,  to  analytically  determine  the 
effects  of  turbulence  on  mobile  parametric  reception,  has  been  accom¬ 
plished.  Theoretical  expressions  have  been  developed  for  the  amplitude 
and  phase  fluctuations  produced  by  turbulence  in  the  interaction  region 
between  pump  transducer  and  hydrophone.  Predictions  have  been  made  for 
the  level  of  fluctuations  that  can  be  expected  in  practical  applications. 
In  this  section,  some  conclusions  are  drawn  from  the  results  of  the 
investigation,  and  future  work  is  discussed. 

The  principal  conclusion  that  can  be  made  is  that  turbulence  can 

produce  significant  variations  in  the  amplitude  and  phase  of  an  acoustic 

signal  detected  by  the  parametric  receiver.  For  a  20  m  long  parametric 

receiver  moving  at  a  speed  of  5  kt,  it  has  been  shown  theoretically  that 

rms  amplitude  variations  exceed  50%  of  the  mean  amplitude,  and  rms  phase 

variations  exceed  0.5  rad,  for  pump  frequencies  of  500  kHz  and  higher. 

It  should  be  noted  that  these  results  are  based  upon  the  assumption  of 

complete  longitudinal  correlation  of  both  amplitude  and  phase  fluctuations 

in  the  region  between  pump  and  hydrophone.  For  array  lengths  that  are 

2 

long  compared  to  the  correlation  distance  of  the  fluctuations,  <B  >  and 

2  ^  ^ 

<SpR>  will  be  less  than  predicted  by  Eqs.  (1)  and  (2). 

A  second  conclusion  that  can  be  drawn  from  the  theoretical  results 
is  that  the  level  of  the  fluctuations  is  highly  dependent  upon  the  inten¬ 
sity  and  geometry  of  the  turbulence  present  in  the  interaction  region. 

2  2 

Specifically,  <B„  >  is  dependent  upon  the  structure  constant  C  [see 
PR  n 

Eq.  (3)],  and  <s|R>  is  dependent  upon  the  mean  square  refractive  index 
variations  <^>  and  the  correlation  distance  a. 
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This  strong  dependence  of  fluctuation  levels  upon  the  turbulence 

parameters  means  that  the  accuracy  of  any  theoretical  prediction  of  <B^  > 

2  PR 

°r  <SPR>  *s  Hmited  by  the  accuracy  of  the  turbulence  parameters  used  in 

making  the  calculation.  These  parameters  are  best  determined  experi¬ 
mentally,  a  point  which  leads  to  the  third  conclusion  resulting  from  the 
investigation,  that  experiments  are  needed  to  determine  the  intensity  and 
geometry  of  turbulence  that  can  be  expected  in  a  practical  application. 

Some  data  regarding  these  turbulence  parameters  should  permit  predictions 
of  the  effects  of  turbulence  on  mobile  parametric  reception  that  would  be 
more  accurate  than  those  presented  in  Section  III. 

Furthermore,  experimental  "testing"  of  the  assumptions  made  in  the 
theoretical  analysis  would  extend  the  usefulness  of  the  present  study  by 
defining  the  limits  of  its  applicability.  For  example,  it  may  be  found 
that  the  assumptions  of  isotropic  turbulence  and  of  complete  transverse 
correlation  of  fluctuations  are  valid  only  for  certain  velocities  or  for 
certain  geometrical  configurations.  As  mentioned  above,  the  longitudinal 
correlation  of  fluctuations  in  the  interaction  region  of  the  parametric 
receiver  may  have  great  impact  upon  the  theoretical  predictions.  Measure¬ 
ments  of  the  longitudinal  correlation  coefficient  for  the  pump  wave  should 
aid  in  the  prediction  of  fluctuation  levels  for  practical  applications. 

Finally,  there  are  two  effects  of  turbulence  that  have  not  been 
considered  in  the  present  study.  One  is  the  spectral  broadening  of  the 
pump  wave,  which  will  contribute  to  the  self-noise  of  the  parametric 
receiver,  and  could  reduce  the  minimum  signal  level  that  it  can  detect. 

A  second  effect  is  the  noise  generated  by  turbulence,  both  volume  dis¬ 
tributed  noise  and  flow  noise  at  the  face  of  the  hydrophone.  This  noise 
could  also  reduce  the  minimum  detectable  level  of  the  parametric  receiver. 
Both  of  these  effects  are  more  amenable  to  experimental  than  to  theoret¬ 
ical  study. 

An  experimental  investigation  generally  dealing  with  the  points 
discussed  above  has  been  proposed  and  will  be  sponsored  by  NAVSEA  Code  63R, 
under  Contract  N00024-79-C-6358,  Task  8. 
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APPENDIX  A 

THE  PARAMETRIC  ACOUSTIC  RECEIVING  ARRAY 
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The  parametric  acoustic  receiving  array  (PARRAY)  is  an  application 

of  the  parametric  array,  which  was  formulated  by  Westervelt^  in  1960  and 

2-12 

has  been  the  subject  of  numerous  studies  in  the  past  two  decades. 

13-19 

Recently  an  integrated  program  at  Applied  Research  Laboratories, 

The  University  of  Texas  at  Austin  (ARL:UT),  has  demonstrated  the  useful¬ 
ness  of  the  PARRAY  as  a  practical  acoustic  sensor. 

Basic  elements  of  the  PARRAY  and  its  operation  are  illustrated  in 
Fig.  A-l.  The  pump  oscillator  and  power  amplifier  generate  the  high  fre¬ 
quency  continuous  signal  that  is  projected  by  the  pump  transducer.  The 
pump  wave  is  shown  in  the  figure  as  closely  spaced,  concentric  arcs. 
Ambient  low  frequency  acoustic  waves,  such  as  the  one  shown  by  widely 
spaced  diagonal  lines,  will  interact  nonlinear ly  with  the  pump  wave  to 
generate  intermodulation  products.  The  function  of  the  receiver  elec¬ 
tronic  hardware  is  to  recover  the  information  contained  in  the  ambient 
signal  by  demodulating  the  interaction  products,  which  appear  as  modula¬ 
tion  sidebands  on  the  pump  carrier.  Thus  an  ambient  acoustic  wave  of 
frequency  f  produces  an  electrical  signal  at  the  output  of  the  receiver 

electronics  which  also  has  frequency  f  . 

s 

The  directional  response  is  identical  to  that  of  a  continuous, 
end-fired  array  of  length  L  equal  to  the  separation  between  pump  and 
hydrophone.  This  directional  response  is  symmetric  about  the  line  joining 
pump  and  hydrophone-  so  that  the  PARRAY  forms  a  conical  beam  with  half 
power  beamwidth  0,  in  degrees,  given  approximately  by 

9  =  105  /X7T 

where  X  is  the  acoustic  wavelength  of  the  signal  to  be  detected. 

A  number  of  desirable  characteristics  result  from  the  fact  that  the 

13 

PARRAY  is  a  continuous,  end-fired  virtual  array.  These  characteristics 
are  summarized  below. 
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FIGURE  A-1 

PARRAY  FUNCTIONAL  DIAGRAM 
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•  Vertical  Directivity  -  Since  the  directional  response  is  symmetric 
about  the  line  joining  the  pump  and  hydrophone,  the  PARRAY  provides  ver¬ 
tical  as  well  as  horizontal  discrimination  against  noise. 

•  No  Grating  Lobes  -  Grating  lobes  are  not  generated  as  the  signal 
frequency  increases  because  the  PARRAY  is  a  continuous  end-fired  array. 

•  Good  Sidelobe  Behavior  -  The  sidelobes  are  well  behaved  and 
decrease  monotonically  to  a  minimum  on  the  back  side  of  the  PARRAY. 

•  High  Front-to-Back  Ratio  -  The  PARRAY  is  relatively  insensitive 
to  signals  arriving  from  the  back  side. 

•  Wide  Bandwidth  -  The  PARRAY  is  inherently  wideband  because  the 
heterodyne  process  translates  the  absolute  bandwidth  of  the  high  fre¬ 
quency  transducers  to  the  low  frequency  signal  region. 

•  Minimum  Number  of  Transducers  -  Two  relatively  small  high  frequency 
transducers  are  required  to  form  the  PARRAY  because  the  nonlinearity  of 
the  water  is  exploited  to  synthesize  the  array  in  the  region  between  the 
transducers . 


APPENDIX  B 

METHODS  FOR  DESCRIBING  TURBULENCE 


Turbulence  has  been  the  subject  of  extensive  study  since  the 

pioneering  work  of  Reynolds  a  century  ago.  Some  excellent  general  ref- 
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erences  on  turbulence  are  the  books  by  Batchelor,  Hinze,  and  Monin 
22 

and  Yaglom.  Introductory  material  may  be  found  in  an  article  by 
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Corrsin  and  in  the  books  by  Tennekes  and  Lumley  and  by  Tritton. 

In  this  appendix,  some  basic  concepts  of  turbulence  are  discussed,  with 

emphasis  on  the  mathematical  methods  used  in  describing  turbulence. 

The  distinguishing  feature  of  turbulent  flow  is  that  the  fluid 
velocity  is  a  random  variable  at  any  given  position  or  time.  The  fluid 
velocity  at  point  r  and  time  t  can  be  represented  as  the  sum  of  mean  and 
fluctuating  components,  i.e., 

v(r,t)  =  v  +  v' 

where 

v  is  the  mean  velocity,  and 

v'  is  the  fluctuation  velocity,  which  is  a  random  variable. 

The  "randomness"  of  the  problem  implies  that  it  will  be  necessary  to 
employ  statistical  parameters  to  describe  the  flow.  One  such  description 
could  be  obtained  by  determining  the  joint  probability  distribution  of 
the  velocity  for  some  number  of  points  in  the  flow.  This  method  becomes 
impractical  if  the  number  of  points  selected  is  large,  and  a  simpler 
description  is  required. 

Although  the  fluid  velocity  is  random,  it  is  not  discontinuous,  so 
there  is  some  correlation  between  the  velocity  at  points  spaced  suf¬ 
ficiently  close  together  in  the  flow.  This  suggests  that  turbulence  can 
be  characterized  by  the  spatial  correlation  coefficient, 


where  the  overbar  denotes  time  averaging.  The  denominator  of  this 

expression  is  just  a  normalizing  factor,  so  R  basically  describes  the 

correlation  of  the  fluctuation  velocity  at  points  r^  and  r^  in  the  flow. 

If  r^55^,  the  correlation  will  be  maximum,  and  R  will  equal  unity.  As 

->■ 

the  separation  between  and  increases,  the  correlation  decreases 

until,  at  very  large  separations,  R  approaches  zero.  A  survey  of  the 

20-25 

literature  related  to  turbulent  flow  shows  that  the  correlation 

coefficient  just  defined  is  one  of  the  most  frequently  used  parameters 
for  describing  turbulence.  With  small  velocity  probes,  such  as  hot-wire 
anemometers,  measurements  of  R  in  turbulent  flow  are  often  made  by  time 

-V  -V 

averaging  the  velocity  fluctuations  at  several  separations  r2~r\' 

There  is  another  method  of  describing  turbulence  that  is  based  upon 
the  kinetic  energy  associated  with  the  flow  rather  than  the  fluid  veloc¬ 
ity.  To  develop  the  concepts  in  an  orderly  fashion,  we  will  need  to  dis¬ 
cuss  the  formation  of  turbulent  eddies. 

One  parameter  associated  with  viscous  fluid  flow  is  the  Reynolds 
number.  Re,  defined  as 


Re 


where 

L  is  the  characteristic  scale  of  flow, 
v  is  the  characteristic  flow  velocity,  and 
v  is  the  kinematic  viscosity  of  the  fluid. 

For  small  values  of  Re,  the  fluid  flow  is  orderly  or  laminar.  When  the 
Reynolds  number  exceeds  a  critical  value,  Recr>  the  flow  becomes  unstable 
and  breaks  up  into  turbulent  eddies.  Each  eddy  of  size  l  will  have  asso¬ 
ciated  with  it  a  local  Reynolds  number,  Re^.  If  Re^  also  exceeds  Recr» 
the  eddy  will  break  down  further  into  smaller  eddies.  This  process  will 
continue  until  the  eddies  are  small  enough  that  viscous  dissipation 
balances  out  the  energy  being  supplied  to  the  eddies  from  the  external 
source.  When  this  occurs,  the  turbulence  reaches  steady  state,  and  a 
range  of  eddy  sizes  exists  such  that 
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Turbulence  has  been  the  subject  of  extensive  study  since  the 

pioneering  work  of  Reynolds  a  century  ago.  Some  excellent  general  ref- 

20  21 

erences  on  turbulence  are  the  books  by  Batchelor,  Hinze,  and  Monin 
22 

and  Yaglotn.  Introductory  material  may  be  found  in  an  article  by 
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Corrsin  and  in  the  books  by  Tennekes  and  Lumley  and  by  Tritton. 

In  this  appendix,  some  basic  concepts  of  turbulence  are  discussed,  with 

emphasis  on  the  mathematical  methods  used  in  describing  turbulence. 

The  distinguishing  feature  of  turbulent  flow  is  that  the  fluid 
velocity  is  a  random  variable  at  any  given  position  or  time.  The  fluid 
velocity  at  point  r  and  time  t  can  be  represented  as  the  sum  of  mean  and 
fluctuating  components,  i.e.. 


v(r , t)  =  v  +  v' 


where 

v  is  the  mean  velocity,  and 

v'  is  the  fluctuation  velocity,  which  is  a  random  variable. 

The  "randomness"  of  the  problem  implies  that  it  will  be  necessary  to 
employ  statistical  parameters  to  describe  the  flow.  One  such  description 
could  be  obtained  by  determining  the  joint  probability  distribution  of 
the  velocity  for  some  number  of  points  in  the  flow.  This  method  becomes 
impractical  if  the  number  of  points  selected  is  large,  and  a  simpler 
description  is  required. 

Although  the  fluid  velocity  is  random,  it  is  not  discontinuous,  so 
there  is  some  correlation  between  the  velocity  at  points  spaced  suf¬ 
ficiently  close  together  in  the  flow.  This  suggests  that  turbulence  can 
be  characterized  by  the  spatial  correlation  coefficient. 


v' (r^)v’ (r2) 
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v'^)  v’(r2) 
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I  <  i  <  L 
o  o 


where 

Lq  is  the  outer  scale  of  turbulence,  and 

£  is  the  inner  scale  of  turbulence, 
o 

The  outer  scale  Lq  is  determined  by  the  boundary  conditions  of  the  fluid 
flow.  In  an  acoustics  application,  Lq  is  generally  taken  to  be  the  di¬ 
mension  from  the  acoustic  source  or  receiver  to  the  nearest  boundary  of 

26  2  7 

the  medium.  The  inner  scale,  i  ,  is  shown  by  Tatarski  ’  to  be 

o 


where  e  is  the  energy  dissipated  as  heat  per  unit  mass  per  unit  time. 

We  see  from  this  expression  that  the  inner  scale  of  turbulence  is  deter¬ 
mined  by  the  viscosity  and  rate  of  energy  dissipation  in  the  medium, 
and  is  independent  of  the  flow  geometry. 

In  describing  turbulence  it  is  often  convenient  to  deal  with  the 
wave  number  k  associated  with  an  eddy  rather  than  its  characteristic 
dimension  l .  The  wave  number  is  inversely  proportional  to  the  eddy  size 
and  may  be  written  as 

2r 

<  =  T  • 

Large  wave  numbers  correspond  to  small  eddy  sizes,  and  vice  versa. 

As  discussed  above,  a  turbulent  flow  with  a  large  Reynolds  number 

will  contain  a  variety  of  eddy  sizes  ranging  from  to  Lq.  Because  the 

kinetic  energy  of  the  flow  is  distributed  throughout  a  spectrum  of  eddy 

sizes,  it  is  possible  to  define  a  power  spectral  density  for  the  flow. 

The  average  kinetic  energy  of  the  flow  is  then  the  sum  of  the  energy 
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associated  with  all  the  eddies,  i.e.,  * 
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a 

T-/ 


4>v(»c)d»c 


where 

T  is  the  average  kinetic  energy  per  unit  mass,  and 

$  (k)  is  the  power  spectral  density. 

^  28 

An  example  of  a  spectrum  for  turbulent  flow  is  shown  in  Fig.  B-l.  In 

the  figure,  the  range  of  wave  numbers  corresponding  to  large,  anisotropic 

eddies  is  labeled  the  "source"  subrange,  as  these  eddies  supply  energy  to 

the  entire  spectrum.  The  energy  of  the  source  eddies  gradually  becomes 

less  anisotropic  (directional)  in  the  "transitional"  subrange  until,  in 

the  "inertial"  subrange,  the  energy  is  isotropic  and  homogeneous.  In  the 

28 

inertial  subrange,  the  spectral  density  is  given  by  the  simple  relation 


0  (k)  =  bx 
v 


-5/3 


9 


where  b  is  a  function  of  the  viscosity  v  and  the  energy  dissipation  rate 
e.  The  -5/3  power  law  in  this  expression  was  originally  proposed  by 
Kolmogorov,  and  has  been  experimentally  verified  in  a  number  of  studies 
(e.g.,  see  Refs.  25-28).  The  turbulent  energy  of  the  flow  is  finally  lost 
to  heat  in  the  "dissipation"  subrange,  where  viscous  forces  become 
dominant . 

In  studying  waves  propagating  in  a  turbulent  medium,  it  is  useful  to 
describe  the  turbulent  field  in  terms  of  its  refractive  index.  When 
turbulence  is  present  in  an  acoustic  medium,  the  sound  speed,  as  a  result 
of  convection  by  the  turbulent  velocity,  will  be  a  random  variable.  The 
sound  speed  may  then  be  written  as 

c  =  c(x,y,z,t)  , 
and  the  refractive  index  is  then 
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n(x,y,z, t) 


c 

o 


c(x,y,z,t) 


1  +  u(x,y,z,t) 


where 

cq  ■  <c>  is  the  mean  sound  speed, 

<  >  denotes  ensemble  averaging,  and 

u(x,y,z,t)  is  the  deviation  from  unity  in  the  refractive  index. 

The  two  methods  of  representing  the  velocity  field  in  turbulent  flow, 
that  of  the  correlation  function  and  that  of  the  spectral  density  func¬ 
tion,  can  also  be  used  to  represent  the  refractive  index  field.  The  cor- 
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relation  function  for  the  refractive  index  variations  is  defined  to  be 


N12  =  u(r2,t) 

If  p(r,t)  is  a  spatially  homogeneous  process,  then  the  correlation  func- 

4-  4- 

tion  depends  only  upon  the  separation  Ar-^-r^ . 

When  the  separation  between  points  becomes  zero,  then  the  correlation 

2 

function  has  a  maximum  of  p  ,  the  mean  square  variation  in  refractive 
index. 

If  the  variations  in  refractive  index  are  isotropic  as  well  as 
homogeneous,  then  the  correlation  function  will  depend  only  upon  the  mag¬ 
nitude  of  the  separation  between  points.  This  condition  may  be  written 
as 


N12  ’  «12M 


where  p=|Ar|.  It  is  useful  to  normalize  the  correlation  function  by 
~ 

dividing  by  P  •  This  results  in  the  correlation  coefficient  R^, 
given  by 
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A  correlaCion  coefficient  frequently  used  to  describe  refractive  index 
variations  is  the  Gaussian  function, 

Ry(p)  =■  exp(-p2/a2)  , 

where  a  is  a  constant  corresponding  to  the  mean  patch  radius. 

As  in  the  case  of  the  turbulent  velocity  field,  the  refractive  index 
variations  can  also  be  described  by  spectral  functions.  These  two  methods 
of  describing  the  random  medium,  that  of  correlation  functions  and  that  of 
spectral  functions,  are  related  by  the  Fourier  transform  theo-em.  In 
three  dimensions,  for  a  homogeneous  medium  this  relation  may  be  written 
as 


N12(Ax,Ay,Az) 


00 

/// 


S  (k., 
y  1’ 


.K3) 


exptj  (K1x+<2y+lc3z^dKld,c2dlc3  ’  t8-1) 


where 

S  is  the  three-dimensional  spectrum  of  the  refractive 

index  variations, 

Kp  <2>  and  are  wave  numbers  in  the  x,  y,  and  z  dimensions, 
respectively,  and 

Ax,  Ay,  and  Az  are  the  x,  y,  and  z  components,  respectively,  of  Ar. 
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For  an  isotropic  medium  Eq.  (B-l)  becomes  « 

00 

N12(p)  =/  V*)dK 

o 
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where  the  one-dimensional  spectral  density  is  related  to  the  Isotropic 
three-dimensional  spectral  density  by 

4>  (k)  =  4iric2S  (ic) 

W  U 

The  spectrum  $  of  the  refractive  index  variations  will  have  the  same 

28 

shape  as  the  velocity  spectrum  4>v,  so  the  characteristics  of  the  sub¬ 
ranges  shown  in  Fig.  B-l  apply  to  both  spectra.  There  is  an  important 

quantity,  the  structure  constant  C  ,  that  can  be  defined  for  the  inertial 

n 

subrange  of  S^.  The  structure  constant  is  a  measure  of  the  intensity  of 
the  refractive  index  variations,  and  its  determination  is  critical  when 
effects  of  turbulence  on  wave  propagation  are  to  be  determined.  In  the 
inertial  subrange,  Cr  is  related  to  the  spectral  function  by 

S  (k)  =  0.333  C2  k“11/3  ,  *  <  k  <  < 

U  n  t  o 

In  Appendix  C,  it  will  be  shown  that  the  parameters  used  to  describe 
turbulence  play  a  large  role  in  determining  the  effects  of  turbulence  on 
a  propagating  acoustic  wave. 


APPENDIX  C 

WAVE  PROPAGATION  IN  INHOMOGENEOUS  MEDIA 


The  effects  of  medium  inhomogeneities  on  the  propagation  of  acoustic 
and  electromagnetic  waves  has  received  a  considerable  amount  of  study  in 
the  past  three  decades.  In  this  appendix  some  of  the  basic  concepts  of 
wave  propagation  in  inhomogeneous  media  are  discussed,  and  pertinent 
theoretical  results  from  the  literature  are  summarized.  Further  details 
may  be  found  in  Refs.  26-33. 

Any  inhomogeneities  in  an  acoustic  medium,  whether  they  be  bubbles, 
biological  matter,  thermal  patches,  or  turbulent  eddies,  may  be  modeled 
as  variations  in  the  sound  velocity  (or  refractive  index)  of  the  medium. 

When  there  are  inhomogeneities  present,  they  scatter  a  propagating  acous¬ 
tic  wave.  Consequently,  the  total  pressure  at  an  observation  point  will 
be  the  sum  of  an  unscattered  pressure  and  a  scattered  pressure.  Since 
the  inhomogeneities  produce  the  scattered  pressure  field,  they  may  be 
treated  as  acoustic  sources.  As  a  first  approximation,  the  inhomogeneities 
are  modeled  as  spherical  sources  of  radius  l,  where  l  is  the  correlation 
distance  of  the  refractive  index  variation.  If  the  dimensions  of  the 
inhomogeneities,  or  "patches,"  are  large  compared  to  the  acoustic  wave¬ 
length,  the  scattered  sound  will  propagate  in  the  same  direction  as  the 

incident  wave  with  a  far field  beam  angle  of  1/kX  radians,  where  k  is  the 

34 

acoustic  wave  number  (see  Fig.  C-l). 

2 

The  region  is  front  of  the  patch  out  to  a  distance  k£  is  the 

nearfield  or  "ray  region"  of  the  scattered  radiation.  In  this  region 

the  patch  will  behave  like  a  lens,  focusing  or  defocusing  the  scattered 

rays  according  to  whether  the  sound  velocity  in  the  patch  is  smaller  or 

34 

larger,  respectively,  than  its  average  value.  Distances  from  the  patch 
2 

farther  than  ki  are  in  the  farfield  or  "wave  region"  of  the  patch. 

One  simple  model  of  an  inhomogenous  medium  is  a  continuous 
distribution  of  patches,  each  patch  having  radius  a  equal  to  the  mean  cor¬ 
relation  distance  of  the  refractive  index  variations.  If  an  acoustic  trans¬ 
mitter  and  receiver  are  located  in  such  a  medium,  the  effect  that  the  medium 

has  on  the  propagating  acoustic  wave  will  depend  upon  the  range  L  of  the 

2  r 

receiver  from  the  transmitter.  For  ranges  less  than  ka  ,  the  receiver  is  in 
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the  nearfield  of  patches  located  in  the  acoustic  path.  In  this  case, 
effects  of  the  inhomogeneities  on  the  propagating  wave  will  be  largely  due 
to  the  phase  delay  and  the  focusing  or  defocusing  that  occurs  as  the  wave 
passes  through  the  patches.  The  range  is  too  short  for  the  path  differences 


between  the  unscattered  signals  and  the  scattered  signals  to  cause  sig¬ 


nificant  diffraction  effects.  Conversely,  for  ranges  Lr>>ka  the  receiver 


is  in  the  farfield  of  most  patches  in  the  acoustic  path.  In  this  case, 
interference  will  occur  between  waves  scattered  by  various  patches,  and 
diffraction  effects  will  dominate.  To  distinguish  between  these  two 
range  conditions,  the  wave  parameter  D  is  used,  and  is  defined  as 


4L 


D  = 


ka 


For  D<<1,  the  receiver  is  in  the  ray  region  or  nearfield  of  the  patches, 
and  for  D>>1,  the  receiver  is  in  the  wave  region  or  farfield. 


The  main  effect  of  medium  inhomogeneities  on  a  propagating  wave  is 
to  produce  fluctuations  in  the  wave’s  amplitude  and  phase.  These  fluc¬ 
tuations  occur  due  to  the  motion  of  the  inhomogeneities,  producing  random 
changes  in  the  scattered  pressure  at  the  receiver.  In  determining  the 
levels  of  fluctuation  in  the  received  wave,  only  a  portion  of  the  in¬ 
homogeneities  in  the  medium  need  to  be  considered.  For  a  plane  wave  (see 

Fig.  C-2(a))  only  the  volume  contained  within  a  cone  of  solid  angle  l/k£, 
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centered  at  the  receiver,  is  significant.  Scattered  pressure  from 

patches  outside  this  cone  does  not  reach  the  receiver.  For  a  spherical 

wave  (see  Fig.  C-2(b))  the  significant  scattering  volume  is  a  narrow 

36 

ellipsoid  with  foci  located  at  the  transmitter  and  receiver.  This 

scattering  volume  is  approximately  contained  within  two  facing  cones  of 

36 

opening  l/k£  centered  at  the  projector  and  hydrophone. 

A  parameter  used  to  characterize  the  amplitude  fluctuations  of  a 
wave  propagating  through  an  inhomogeneous  medium  is  the  coefficient  of 
amplitude  variation,  CAV,  CAV  is  defined  as  the  standard  deviation  in 
amplitude.  A,  of  the  wave,  and  is  written  as 


HYDROPHONE 


(a)  FOR  PLANE  WAVE 


(b)  FOR  SPHERICAL  WAVE 


FIGURE  C-2 
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CAV  = 


2  2  1/2 
CA2-Ao2) 


where  A  =A  and 
°35 

reasons,  the 
defined  as 


the  overbar  indicates  time  averaging.  For  mathematical 
log  amplitude  fluctuation  B  is  often  used,  where  B  is 


B  = 


For  small  fluctuations, 

A  -  A 

B  =  a"£  ,  B  «  1 

o 

~~2  2 

and  B  is  equivalent  to  CAV  .  The  phase  fluctuations  of  the  wave  are 
represented  by  the  mean  square  phase  deviation,  where  the  phase  deviation 
is 


S  =  $  -  4>0 

<J>  is  the  instantaneous  phase,  and  <j>  =  (j>.  The  mean  square  phase  deviation 

o 

is  therefore 


Often  the  ensemble  averages,  denoted  by  <B  >  and  <S  >,  are  used 

~2  ~~2 

rather  than  the  time  averages  B  and  S  .  For  ergodic  processes  these  two 
types  of  averages  are  equivalent. 


Theoretical  results  have  been  obtained  for  the  amplitude  and  phase 
fluctuations  of  waves  propagating  through  a  statistically  isotropic 
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1 


inhomogeneous  medium.  For  large  values  of  the  wave  parameter  (D»l) , 
35 

Chernov  obtains 


<B2> 


2 

<S  > 


2 

<P  > 


k2L 


00 


R(p)dp 


(C-l) 


where 

p  is  the  deviation  from  unity  in  the  refractive  index  n,  i.e., 
n  =  c q/c  -  1  +  p, 

cq  is  the  average  of  the  sound  speed,  c,  and 
R(p)  is  the  correlation  coefficient  of  the  refractive  index 
variations. 

2  2 

If  the  correlation  coefficient  is  Gaussian,  given  by  exp(-p  /a  ),  then 
Eq.  (C-l)  becomes 


<B2>  =  <S2>  =  y  <p2>  k2aLr  •  (C-2) 
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Mintzer  has  shown  ’  that  Eqs.  (C-l)  and  (C-2)  are  valid  for  a  spherical 
wave  in  the  wave  region  when 


2  2 

k  <  p  >  aL  <<  1 

r 

For  small  values  of  the  wave  parameter  (D<<1)  corresponding  to  the 

ray  region  of  the  patches,  Chernov  shows  the  mean  square  amplitude  fluc- 

35 

tuations  for  a  plane  wave  to  be 


<B2> 


1  2 
=  -  <p  > 


C( 

if 


2  2 
V  V 


R(p)dp 


D«1 


which,  for  a  Gaussian  correlation  coefficient,  becomes 


(C-3) 
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,  D«1 


(C-4) 


<b2> 


8/ir  2 

—  <P  > 


_r 

3 
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The  phase  fluctuations  for  D<<1  have  been  shown  to  be  twice  their 
values  when  D>>1  [i.e.,  double  the  results  in  Eqs.  (C-l)  and  (C-2) ] . 


Finally,  when  the  wave  parameter  is  of  the  order  of  unity,  Chernov 

L  ,  35 

obtains 


_2  /n  2  ,  2  T 

<B  >  =  —  <y  >  k  a  L 


t(1  -5 


—  arctanD 


/ 


(C— 5 ) 


<S2>  =  <y2>  k2a  Lr|l  +  -jj  arctanDj  , 


(C-6) 


assuming  a  Gaussian  correlation  coefficient.  For  small  and  large  values 

of  D,  these  last  two  results  become  equivalent  to  the  ray  and  wave  region 

solutions,  respectively.  It  is  assumed  in  obtaining  the  above  results 

that  the  inhomogeneities  are  large  compared  to  an  acoustic  wavelength 

(ka>>l),  and  that  the  propagation  distance  is  large  compared  to  the  scale 

of  the  inhomogeneities  (L  >>a) .  Using  the  same  conditions  and  assump- 
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tions  as  Chernov,  Karavainikov  derived  a  similar  set  of  results  for  the 
amplitude  and  phase  fluctuations  in  a  propagating  spherical  wave.  For 
the  wave  region  of  the  inhoraogeneities  where  D>>1,  Karavainikov 's  result 
is  identical  to  Eq.  (C-2). 


The  results  discussed  above  pertain  to  a  medium  that  can  be  modeled 
as  containing  a  continuous  distribution  of  patches  of  radius  a.  As  dis¬ 
cussed  in  Appendix  B,  a  turbulent  medium  contains  a  distribution  of 
eddies  ranging  in  size  from  the  inner  scale  £  to  the  outer  scale  Lq. 

For  the  case  of  a  medium  whose  refractive  index  field  is  determined  by 

26  27  31 

turbulence,  techniques  have  been  developed  by  Tatarski  ’  and  others 

for  calculating  wave  parameter  fluctuations  similar  to  the  results  sum- 

26  27 

marized  above  for  a  Gaussian  medium.  In  particular,  Tatarski  *  has 
found  expressions  for  amplitude  fluctuations  in  terms  of  a  turbulence 
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parameter,  C^.  For  propagation  distances  such  that  Lo>>/ALr>>!0, 
where  A  is  the  acoustic  wavelength,  Tatarski  shows 


<B2>  ■»  0.13  C2  k 7/6  L11/6 
n  r 


(C-7) 


for  spherical  waves,  and 


0.31  C2  k7/6  LU/6 
n  r 


(C-8) 


for  plane  waves.  The  quantity  Cn  is  related  to  the  spectral  function 
S^(k)  (defined  in  Appendix  B)  by^ 

S  («)  =  0.033  C2  kt-11/3  ,  x  <  <  <  k 

u  n  t  p 

An  assumption  made  in  deriving  Eqs.  (C-7)  and  (C-8)  is  that  only  the 
inertial  range  of  the  spectrum  contributes  to  the  acoustic  fluctuations. 
If  the  spectrum  is  of  the  form  shown  in  Fig.  B-l,  so  that  there  is  an 
approximately  flat  transition  region,  then  there  is  an  additional  contri¬ 
bution  from  the  transition  region  to  the  amplitude  fluctuations  of  an 
amount 


<B2>  =  -r^77  $  L 
t  480  m 


(C-9) 


for  spherical  waves,  and 


<B2>  =  jz  <t>  L3(<^-x4) 

t  48  m  r\  t  m/ 


(C-10) 
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for  plane  waves.  As  shown  in  Fig.  B-l,  4>m  is  the  maximum  value  of  the 
spectral  density  function  4>^.  The  total  mean  square  amplitude  fluctua¬ 
tion,  when  there  is  a  contribution  due  to  the  transition  region,  is 

2  2 
<B  >+<b£> . 


42 


J 


The  results  summarized  In  this  appendix  are  used  In  Appendix  D  to 
obtain  expressions  for  the  amplitude  and  phase  fluctuations  in  the  side¬ 
band  pressure  of  the  parametric  receiver. 


43 


THEORETICAL  ANALYSIS 
IN  A 


Itffjtf;) 


A. 


Introduction 


The  presence  of  turbulent  eddies  in  the  interaction  region  of  a 
mobile  parametric  receiver  will  have  adverse  effects  on  receiver  perform¬ 
ance.  Some  effects,  specifically  the  fluctuations  in  the  amplitude  and 
phase  of  the  detected  acoustic  signal,  are  analyzed  in  this  appendix. 


The  approach  taken  in  this  analysis  is  to  treat  the  turbulent  eddies 
as  inhomogeneities  that  scatter  acoustic  waves  propagating  through  the 
turbulence,  as  shown  in  Fig.  D-l.  When  inhomogeneities  are  present  in 
the  interaction  region,  the  virtual  sources,  qdv,  depend  not  only  upon 


the  direct  radiation  p  and  p  from  the  pump  and  signal  sources, 

p  s 


respectively,  but  also  upon  the  scattered  pump  pressure,  p^.  As  the 

scatterers  move  about  in  time,  the  phase  difference  between  p  and  p' 

P  P 

changes,  causing  amplitude  and  phase  variations  in  the  source  density 


function,  q.  Similarly,  the  radiation  from  the  virtual  sources  is 
scattered,  causing  further  fluctuations  in  the  amplitude  and  phase  of  the 
sideband  signals  detected  by  the  hydrophone. 


B.  The  Second-Order  Pressure  at  the  Hydrophone 

The  first  major  step  in  the  analysis  is  to  develop  an  expression  for 
the  second-order  pressure  field  produced  by  nonlinear  interaction  of  the 
pump  and  signal  waves.  This  can  be  done  by  accounting  for  the  fluctua¬ 
tions  in  the  pump  and  signal  wave  that  are  produced  by  turbulence,  and 
then  calculating  the  Westervilt  source  density  function. 


The  variations  in  the  amplitude  and  phase  of  the  pump  wave  are  taken 
into  account  by  writing 


„  -a  r  j  (u»  t-k  r-S  ) 

p  =  j(P  /r)  D(y)(l+B  )  e  P  e  P  P  p  ,  (D-l) 

P  P  P 
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where 


Pp  is  the  mean  pump  pressure  amplitude  at  Y*0,  r=l  m, 

D(y)  is  the  directivity  function  for  the  pump  transducer, 


B  =  [P  (t)-P  ]/P  is  the  fractional  variation  of  the  instantaneous 
p  P  P  P 

pressure  amplitude  P  (t)  from  its  mean  value,  and 


4> (t )  — (t)  is  the  variation  of  the  instaneous  phase  $(t)  of  the 


pump  wave  from  its  mean  value,  $(t). 

The  remainder  of  the  expression  represents  a  spherical  wave  propagating 
radially  outward  from  an  origin  situated  at  the  pump  (see  Fig.  D-2),  with 
angular  frequency  to^,  wave  number  k^,  and  attenuation  coefficient  a^. 


Inclusion  of  the  attenuation  coefficient  as  a  constant  requires  that 
the  time  variations  in  amplitude  and  phase  are  slow  compared  to  the  pump 
frequency.  In  other  words,  frequency  broadening  of  the  pump  wave  due  to 
the  fluctuations  is  small  enough  that  the  attenuation  coefficient  remains 
constant. 


Scattering  of  the  signal  wave  in  the  interaction  region  is  assumed 

to  be  negligible  compared  to  scattering  of  the  pump  wave  because  of  the 

difference  in  frequency  of  the  two  waves.  The  signal  wave  is  therefore 

modeled  as  a  plane  wave  having  frequency  u>  wave  number  k  ,  and  attenua- 

S  y  s 

tion  coefficient  a  : 

s 


P 


s 


P 

s 


-a  z 
s 
e 


j(“st_ksz) 

e 


(D-2) 


where  Pg  is  the  pressure  amplitude,  a  constant.  It  has  been  assumed  for 
convenience  that  the  acoustic  source  generating  pg  is  located  on  the 
maximum  response  axis  of  the  parametric  receiver,  as  shown  in  Fig.  D-2. 
The  signal  source  is  assumed  to  be  sufficiently  far  from  the  pump  trans¬ 
ducer  that  the  signal  wave  is  approximately  planar  in  the  vicinity  of  the 
parametric  receiver. 

The  source  density  for  the  second-order  radiation  can  be  found  using 
the  equation 
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uO>  > 


q(r,t)  = 


(D-3) 


6  3  r  ,  2, 

2  4  3t  1 

P  c 
o  o 


where 

r  is  the  position  vector  of  a  source  point,  as  shown  in  Fig.  D-2, 

6  is  a  parameter  of  nonlinearity, 

Po  is  the  ambient  density, 

cq  is  the  ambient  sound  speed,  and 

p  =  p  +  p  is  the  total  first-order  pressure  field. 

X  p  s 

Substitution  of  Eqs.  (D-l)  and  (D-2)  for  the  total  first-order  pressure 
in  Eq.  (D-3)  gives  the  source  density  function  as 


-6P  P  a)  -a  z  -a  r  -jS  -j  (k  r±k  z)  jw.t  . . 

■  ,-aJ-  «  S  «  P  0<Y)  (l+»  )  .  p  e  P  8  «  1  .  <D-‘> 

2  4  P 

per  r 

o  o 


where 


(j.  =  w  ±  u  ,  and 

ri  s 

r  =  r  . 


It  is  assumed  that  the  pump  wave  fluctuation  terms  B  and  S  have 

P  P 

complete  spatial  correlation  along  the  spherical  wave  fronts  within  the 

pump  beam,  so  that  and  are  independent  of  y.  With  this  assumption, 

a  solution  for  the  second-order  pressure  at  the  hydrophone  can  be  obtained 

9 

by  adopting  a  procedure  developed  by  Berktay  and  Shooter.  They  assume 
that  the  sphericity  of  the  pump  wave  within  the  beam  is  small  compared  to 
the  wavelength  at  the  signal  frequency  so  that,  when  the  signal  source  is 
collinear  with  the  pump  and  hydrophone,  the  virtual  sources  q+dv  can  be 
assumed  to  be  cophasal  on  the  spherical  wave  fronts.  The  frequency  of 
the  second-order  radiation  is  nearly  equal  to  the  pump  frequency,  so  the 
second-order  waves  will  radiate  spherically,  with  the  same  beam  pattern 
as  for  the  pump  wave.  These  assumptions  are  used  to  calculate  the  second- 
order  pressure  at  the  observer  as  follows. 
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The  elemental  particle  velocity  at  r  due  to  a  spherical  shell  of 
sources  of  thickness  fir  Is 


fiur  “  j  q(r)fir 


where  the  source  density  function  is  given  by  Eq.  (D-4).  The  contribu¬ 
tion  of  these  sources  to  the  particle  velocity  at  (L ,0)  will  be 


fii^  *  (r/L) fiu^ (1+B^)  e  U  exp [-(a++jk+) (L-r) ]  , 


where 


Bu  ■  ^  ^ ^  is  the  amplitude  fluctuation, 

S  =  d>  (t)  -  k  r  is  the  phase  fluctuation, 
u  u  ± 

U(t)  and  are  the  amplitude  and  phase  of  the  particle  velocity, 

respectively, 

k+  is  the  acoustic  wave  number  at  frequency  u>+,  and 

a+  is  the  attenuation  coefficient  at  frequency  w+. 

It  is  assumed  that  frequency  broadening  of  the  interaction  frequency  wave 
due  to  the  inhomogeneities  is  small  enough  that  a+  remains  constant.  The 
total  particle  velocity  at  (L,0)  is  the  sum  of  contributions  from  all 
sources  in  the  interaction  region,  i.e.. 


L 

U±(L,0)  =  (2L)_1  exp[-(a±+jk±)L] J (1+Bu) 


e  U  q (r) r  exp[ (a++jk+)r]dr  .  (D-5) 


The  second-order  pressure  at  the  point  (L,0)  can  now  be  found  by  using 
the  farfield  relation,  p=pQcu,  in  connection  with  Eq.  (D-5).  The  result 
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-gP  P  w+ 

p+(L,0)  =  - S  exp[-(ot++jk+)L] 

2p  c  L 
o  o 


L 

/ 


(l+B  )  (1+B±)  e 


'1(VS±) 


X 


exp[-(a  +a  -a  )z]dz 
s  p  t 


(D-6) 


where  B+  and  S+  are  amplitude  and  phase  fluctuations,  respectively,  for 
the  interaction  frequency  pressure  wave.  (It  may  be  shown  that  the 
amplitude  fluctuations  in  the  pressure  and  particle  velocity  are  equal; 
i.e.,  Bu=B+.  Similarly,  Su=S+.)  Since  the  second-order  radiation  is 
cophasal  along  the  spherical  wave  fronts,  a  change  of  variable  from  r  to 
z  has  been  made  in  Eq.  (D-6).  The  physical  interpretation  of  this  result 
follows.  Nonlinear  interaction  of  the  pump  and  signal  waves  (p^  and  pg) 
will  produce  at  the  hydrophone  of  the  parametric  receiver  a  pressure  wave 
p+.  The  amplitude  of  p+  increases  linearly  with  the  parameter  of  non¬ 
linearity,  the  signal  and  pump  pressure  amplitudes,  and  the  sideband  fre¬ 
quency.  Although  the  array  length  L  appears  in  the  denominator  of 
Eq.  (D-6),  some  care  is  needed  in  interpreting  the  dependence  of  p+  upon 
L.  In  a  homogeneous  medium  the  fluctuation  terms  become  zero,  and  since 


the  integral  has  a  value  L  which  cancels  with  the  L  in  the  denominator  of 
Eq.  (D-6).  Thus,  in  the  homogeneous  case  the  amplitude  of  p+  is  inde¬ 
pendent  of  array  length.  When  inhomogeneities  are  present,  both  amplitude 
and  phase  of  p+  depend  upon  fluctuations  caused  by  turbulent  scattering 
in  the  interaction  region  between  pump  anJ  hydrophone.  Ignoring  attenua¬ 
tion,  the  integral  in  Eq.  (D-6)  represents  the  sum  of  the  amplitude  and 
phase  fluctuations  associated  with  each  elemental  length  dz  of  the  para¬ 
metric  receiver. 
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C.  Fluctuations  in  the  Second-Order  Pressure 

Now  that  an  expression  has  been  obtained  for  the  second-order 
pressure  p+  in  a  turbulent  medium,  the  next  step  of  the  analysis  is  to 
find  the  level  of  amplitude  and  phase  fluctuations  in  this  second-order 
pressure.  This  can  be  accomplished  as  follows. 

In  addition  to  Eq.  (D-6) ,  the  sideband  pressure  can  be  written 
explicitly  in  terms  of  the  amplitude  and  phase  fluctuations  observed  at 
the  hydrophone.  If  p^  is  the  sideband  pressure  that  would  be  observed 
in  a  homogeneous  medium  (i.e.,  with  no  turbulence),  then  the  sideband 
pressure  in  the  presence  of  inhomogeneities  is 

~jSPR 

p±  =  Ph(1+BPR)  6  *  (D~7) 

where  B  and  S  account  for  the  total  amplitude  and  phase  fluctuations 
PR  PR 

in  the  sideband  pressure  at  the  hydrophone. 

A  solution  for  the  amplitude  and  phase  fluctuations  B  and  S  can 

PK  PK 

be  obtained  as  follows.  As  a  first  approximation,  it  is  assumed  that  the 
acoustic  wave  fluctuations  are  sufficiently  small  that 


B 

«  1, 

S 

«  1 

P 

P 

B+ 

«  1, 

S+ 

«  1. 

Then  the  amplitude  fluctuation  terms  can  be  expanded  and,  retaining  only 
the  first-order  fluctuations  terms, 

(1+B  ) (1+B  )  =  1  +  B  +  B ,  .  (D-8) 

P  ±  P  ± 

Similarly,  the  phase  fluctuations  can  be  approximated  as 


CD— 9) 


-JCSp+V 


i  1  - 


j (Sp+S±) 


Using  Eqs.  (D-8)  and  (D-9)  for  the  fluctuation  terms,  Eqs.  (D-6)  and 
(D—  7 )  may  be  set  equal,  giving 


Ph(1+BPR) 


-jS 


PR 


=  A 


JU 

/ 


(1+B  +B±)[1  - 


j (sp+s+) ldz 


(D-10) 


where 


-6PsV±  -(a++jk±)L 
A  = - ^ —  e 

2p  c  L 
o  o 

and  the  attenuation  terms  in  the  integrand  of  Eq.  (D-6)  have  been  approxi¬ 
mated  as 


If  the  fluctuations  in  the  sideband  pressure  are  also  assumed  to  be  small, 
then 


a«PR)  e  JSpR^  a+Bp^a-jSpR)  . 

Using  this  approximation  and  retaining  fluctuation  terms  to  first  order 
only,  Eq.  (D-10)  becomes 


Ph(1+BpR_jSPR) 


=  A 


L 

J (l+B-jS)dz 
o 


(D-ll) 
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where 


B  ■  B  +  B  ,  and 
P  ± 

S  =  S  +  S^. 

P  ± 

Subtraction  of 


L 


o 


from  both  sides  of  Eq.  (D-ll)  leaves 


BPR  "  ^  SPR 


S)dz 


Equating  real  and  imaginary  parts  of  this  result  gives 


PR 


x. 

’if 


B  dz 


and 


L 

■if 


SpR  =  f  I  S  dz 


As  B  and  S0  are  random  quantities  with  zero  mean,  it  is  useful  to  deal 
PR  PR 

with  the  mean  squared  amplitude  and  phase  fluctuations, 


L  L 


<BDO> 

PR 


=  J  f  /<BlB2>  dZldZ 5 


(D-12) 


o  o 
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and 


<S1S2>  dz1dz2  ,  (D-13) 

o  o 

where  the  angular  brackets  <>  denote  ensemble  averaging.  It  is  easier 
to  interpret  these  results  physically  by  writing  Eqs.  (D-12)  and  (D-13) 
explicitly  in  terms  of  the  pump  and  sideband  fluctuations.  For  example, 
Eq.  (D-12)  can  be  written  as 


<SPR>  = 


L  L 

?// 


L  L 


<B2  >  = 
PR 


?// 

o  o 


("W  +  <BplB±2> 


+  <B±1B  .>  +  <B+1B±2>)dz1dz2  ,  (D-14) 


where 

the  subscript  p  denotes  the  pump  wave, 

the  subscript  ±  denotes  the  interaction  (or  sideband)  frequency  wave 
radiated  by  the  virtual  sources,  and 
the  subscripts  1  and  2  refer  to  virtual  sources  located  at  z^  and 

Z2- 

The  terms  in  the  integrand  of  Eq.  (D-14)  describe  the  spatial  correlation 
of  the  fluctuations.  For  example,  <BplBp2>  gives  the  correlation  between 
the  amplitude  fluctuations  in  the  pump  wave  at  z^  and  the  amplitude  fluc¬ 
tuations  in  the  interaction  frequency  wave  at  z^.  The  equation  thus 
relates  the  amplitude  fluctuations  in  p+  at  the  hydrophone  to  the  cor¬ 
relations  of  fluctuations  in  pp  and  p+  in  the  interaction  region. 

D.  Approximations  for  the  Spatial  Correlation  Terms 

Evaluation  of  Eq.  (D-14)  requires  that  the  correlation  terms  in  the 
integrand  be  expressed  as  functions  of  z^  and  z2»  These  correlation 
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terms  fall  into  two  categories:  (1)  crosscorrelation  between  pump 
amplitude  fluctuations  and  interaction  frequency  amplitude  fluctuations, 
and  (2)  autocorrelation  of  pump  and  interaction  frequency  amplitude  fluc¬ 
tuations.  These  categories  are  discussed  separately  below. 


(1)  Crosscorrelation.  The  term  <Bp(z1)B+(L-z2)>  represents  the 
crosscorrelation  between  pump  pressure  amplitude  fluctuations  at  the 
point  (z^)  and  interaction  frequency  amplitude  fluctuations  at  (L-z2) . 
The  frequencies  of  the  two  waves  are  approximately  equal  (Wp=u)+)  so 
complete  frequency  correlation  may  be  assumed  between  B^  and  B+. 


The  propagation  paths  associated  with  the  pump  and  interaction 
frequency  waves  are  shown  in  Fig.  D-3.  The  term  B^(z^)  is  due  to  scat¬ 
tering  of  the  pump  wave  as  it  propagates  from  the  origin  to  z=z^.  The 
fluctuations  in  the  interaction  frequency  wave  B+(L-z2)  are  due  to 
scattering  of  the  second-order  radiation  as  it  propagates  from  a  source 
at  point  z2  to  the  hydrophone  at  z=L.  Both  z^  and  z2  may  vary  between 
0  and  L,  so  there  will  be  situations  in  which  the  two  propagation  paths 
overlap  (z^>z2)  and  situations  in  which  they  are  separate  (z^<z2) . 


An  exact  solution  to  Eq.  (D-14)  requires  an  expression  for  the 
correlation  term  <B  B  >  for  all  values  of  z,  and  z„,  and  for  both  near- 

pi  l  4 

field  and  farfield  receiving  arrays.  This  is  essentially  a  problem  of 
calculating  the  correlation  of  amplitude  fluctuations  at  two  receivers 
when  there  are  two  sources  generating  separate  waves.  The  solution  to 
this  problem  is  not  available  in  the  literature,  nor  is  it  readily 
obtained,  but  it  is  convenient  to  make  the  following  simplified 
approximation. 


Two  asymptotic  conditions  are  considered,  one  in  which  the  array 
length  L  is  much  less  than  the  longitudinal  correlation  distance  1^  and 
the  other  in  which  the  array  length  is  much  greater  than  the  pump  cor¬ 
relation  distance.  First  assume  that  L<<1  .  In  this  case  the  separation 
between  the  "receivers"  at  z  and  L  is  less  than  the  pump  correlation 
distance,  so 
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(D-15) 


L  -  z,  <<  1 
1  P 

where  1  is  the  distance  at  which  the  correlation  coefficient  for  the 
P 

pump  wave  amplitude  fluctuations  equals  1/e.  When  Eq .  (D-15)  applies, 
the  amplitude  fluctuations  of  the  pump  wave  at  z^  and  the  interaction 
frequency  wave  at  z=L  will  be  highly  correlated.  Their  correlation 
coefficient  may  be  approximated  as  unity;  i.e.. 


<Bp(z1)B±(L-z2)> 

2  2  f/2 

<Bp(z1)><B+(L-z2)>| 


L  «  1 

P 


From  this  expression  the  crosscorrelation  between  pump  and  interaction 
frequency  fluctuations  can  be  written  as 


<Bp(Zl)B±(L-z2)>  =  |<Bp(z1)xB2(L-z2)>j 


1/2 


(D-16) 


This  equation  is  in  terms  of  mean  square  amplitude  fluctuations,  and  can 
be  calculated  using  the  results  discussed  in  Appendix  C.  A  similar 
approximation  can  be  made  for  the  remaining  crosscorrelation  term  in 
Eq.  (D-14) ,  namely 


<Bp(z2)B±(L-z1) 


(D-17) 


It  should  be  noted  that  Eqs.  (D-16)  and  (D-17)  apply  only  for  short  array 

lengths  (L<<1  ).  In  general,  if  the  array  length  is  such  that  L<<1 
P  p 

then  Rp  +  is  less  than  unity,  and  the  following  less  restrictive  approxi¬ 
mations  apply: 
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(D-18) 


<Bp(z2)B+(L-z1)>  =  [ <  Bp (z2) ><B^(L-z1) >| 


1/2 


and 


1/2 


2  2 

<Bp(z1)B±(L-z2)>  =  <Bp(z1)><B±(L-z2)>  Rp  +  .  (D-19) 


Now  consider  the  situation  when  the  array  length  is  much  greater 
than  the  pump  longitudinal  correlation  distance  (L>>lp) .  Most  separa¬ 
tions  L-z^  will  be  greater  than  the  correlation  distance  1  ;  therefore 
the  fluctuations  will  (on  the  average)  have  very  little  correlation. 
For  L*>lp  the  correlation  coefficient  may  therefore  be  approximated  by 
zero : 


<Bp(Zl)B±(L-z2)>  =  p,  L»1 


(D-20) 


and 


<Bp(z2)B±(L-Zl)>  *  °»  L>>1p  ‘  (D-21) 

These  results  are  used  in  Eq.  (D-14)  for  long  array  lengths. 

(2)  Autocorrelation.  The  second  category  of  correlation  terms  in 
Eq.  (D-14)  are  autocorrelation  functions  for  the  pump  and  interaction 
frequency  amplitude  fluctuations. 


The  pump  autocorrelation  function  will  be  given  by 


<B‘ (z.)B(z,)> 

p  1  p  2 


(D-22) 


where  Rp  is  the  longitudinal  correlation  coefficient  for  pump  wave  fluc¬ 
tuations  at  z^  and  z^.  For  short  lengths  (L<<1  ) ,  Rp  can  be  approximated 
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as  unity.  For  longer  array  lengths  the  results  of  Chernov  or 
42 

Eliseevnin  can  be  used  to  estimate  R  . 

P 

The  autocorrelation  function  f^r  the  interaction  frequency  wave  is 
different  from  that  for  the  pump  wave.  Rather  than  originating  at  a  com¬ 
mon  source  point  and  being  received  at  different  observation  points,  as 
is  the  pump  wave,  the  interaction  frequency  waves  originate  at  different 
source  points  and  are  received  at  a  common  observation  point  as  shown  in 
Fig.  D-4.  The  waves  originate  at  z^  and  z 2  and  are  both  received  at  z=L. 

While  there  is  no  explicit  analysis  of  this  situation  in  the  literature, 

40 

Chotiros  and  Smith  have  demonstrated  that  the  principle  of  reciprocity 
applies  as  follows.  If  a  wave  of  frequency  w+  is  projected  from  the 
transducer  at  z=L,  then  the  correlation  of  amplitude  fluctuations  received 
at  points  z^  and  z 2  will  be  <Br (L-z^) B^ (L-z^) > .  By  the  reciprocity  prin¬ 
ciple,  this  correlation  will  be  identical  to  that  for  waves  originating  at 
z^  and  z^  and  received  at  z=L;  i.e., 

<B±(L-z1)B+(L-z2)>  =  <Br(L-z1)Br(L-z2)> 


where  the  subscript  r  indicates  that  the  positions  of  sources  and  receivers 
have  been  interchanged.  This  result  is  useful  because  it  leads  to 


r  2  2  il/2 

<B±(L-z1)B±(L-z2)>  =  |<B+(L-z1)><B±(L-z2)>]  Rr  ,  (D-23) 


where  R^  is  the  longitudinal  correlation  coefficient  for  interaction  fre¬ 
quency  waves  originating  at  z=L  and  received  at  points  z^  and  z^.  When 

the  array  length  is  short  (L<<1  ),  the  correlation  coefficient  is  R  =1. 

P  41  42  r 

For  longer  array  lengths,  the  results  of  Chernov  or  Eliseevnin  may  be 

used  to  calculate  R^.  Because  complete  frequency  correlation  has  been 

assumed  for  the  pump  and  interaction  frequency  waves,  the  coefficient  R^ 

will  be  equal  to  the  longitudinal  correlation  coefficient  for  a  wave  at 

frequency  on  received  at  ranges  L-z^  and  L-z^i  that  is, 

Rr  =  Rp(L-z1,  L-z2) 
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FIGURE  D-4 

GEOMETRY  FOR  AUTOCORRELATION  OF 
SECOND-ORDER  WAVE  FLUCTUATIONS 


ARL:UT 
AS- 80-1005 
CRC-GA 
3-31-80 


The  mean  squared  fluctuation  terms  in  Eq.  (D-23)  are  for  waves  of 
frequency  u>+  traversing  paths  of  lengths  L-z^  and  L-z2>  and  may  be  cal¬ 
culated  using  the  formulae  discussed  in  Appendix  C. 

Having  obtained  expressions  approximating  the  various  correlation 
terms,  it  is  now  possible  to  determine  the  level  of  amplitude  fluctua¬ 
tions  produced  by  turbulence. 

E.  Amplitude  Fluctuations 

The  mean  squared  amplitude  fluctuations  in  the  second-order  pressure 

2 

detected  by  the  hydrophone,  <BpR>,  can  now  be  written  as 


L  L, 


2  l  r  r  1 1/2  11/2 

<BpR>‘?  JJ  V  <Bp<*l>><Bi<L-z2»  Rp, 

o  O  j  L 


<Bp(z2)><B^(L-Zl)> 


1/2 


R  .  + 
P>± 


2  2  ]1/2 
<B^(L-z1)><B^(L-z2)>  Rx 


(D-24) 


dzrdz2 


where  Eqs.  (D-18),  (D-19) ,  (D-22) ,  and  (D-23)  have  been  substituted  into 
Eq.  (D-14) . 


As  an  example  of  evaluating  Eq.  (D-24),  assume  the  array  length  is 
sufficiently  short  that 


If  it  is  further  assumed  that 

L  >>  A  L  »  S, 

o  p  o 


(D-25) 


(D-26) 


so  that  the  inertial  subrange  is  dominating  the  scattering  effects,  then, 
from  Eq.  (C-7) , 
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( 


<B2>  =  0.13  c2  k7/6  z11/6 
P  n  p 


<B2>  *  0.13  C2  k7/6  (L-z)11/6 
±  n  ± 


Substituting  Eqs.  (D-27)  and  (D-28)  into  (D-24)  gives 


n  n  r2  t7/6  L 

V  .  “  n  %  f f f  11/12  11/12  ,  11/12  .T  .11/12 

<BPR> - p  y  y  Z1  z2  +  Z1  (L~z2) 


(D-27) 


(D-28) 


.  ,11/12  .  .11/12  ..  .11/12  /T  .11/12  .  , 

+  z2  (L-z^^)  +  (L-z^)  (L-z2)  dz1,dz2 


where  it  has  been  assumed  that  k^=k+.  Integration  is  straightforward, 
and  the  mean  squared  amplitude  fluctuations  in  the  second-order  pressure 


<B2  >  =  0.1415  C2  k7/6  L11/6 
PR  n  p 


(D-29) 


This  equation  is  one  of  the  principal  results  of  this  analysis,  for  it 
represents  the  mean  squared  amplitude  fluctuations  in  the  "output"  (upper 
or  lower  sideband  pressure)  of  a  mobile  parametric  receiver  when  turbu¬ 
lence  is  present  between  pump  and  hydrophone.  The  limitations  imposed  on 
this  result  by  the  assumptions  used  in  the  analysis  are  discussed  below. 

F.  Phase  Fluctuations 


An  expression  for  phase  fluctuations  similar  to  Eq.  (D-29)  can  be 
obtained  as  follows.  The  integral  expression  for  the  mean  squared  phase 
fluctuations,  Eq.  (D-13) ,  may  be  expanded  to  give 


>  i 


<SPR>  “  -J  JJ  (<SplSp2>  +  <SplS±2>  +  <S±1SP2>  +  <S±lS±2>jd‘ldz2  •  ( 


The  spatial  correlation  terms  in  the  integrand  of  Eq.  (D-30)  may  be 
expressed  in  terms  of  parameters  characterizing  the  turbulent  medium. 
As  discussed  in  Appendix  C,  the  mean  squared  phase  fluctuations  for  a 
wave  propagating  a  distance  Lr  in  an  inhomogeneous  medium  is  given  by 


<S2>  =  <y2>  k2  Lr  J  R(p)dp 


(C-l) 


where  it  is  assumed  that  the  wave  parameter  is  much  greater  than  unity 

(D>>1) .  Equation  (C-l)  can  be  used  to  obtain  approximate  expressions 

for  the  spatial  correlation  terms  appearing  in  Eq.  (D-30).  For  example, 

if  the  term  <S  ,S  „>  is  written  in  the  form 
pl  p2 

,  1/2  2  1/2 

<VV  ‘  ‘V  ‘V  Ri2  • 


where  R  is  the  longitudinal  correlation  coefficient  for  the  phase  flue 
11  2  2 
tuations,  then  Eq.  (C-l)  can  be  used  to  rewrite  <S  >  and  <Sp2>s  giving 


<S  ,S  0>  = 
pl  P2 


2  .  2 
<p  >  k 

P 


00 

I  R(p)dp  z* 

n  J 


1/2  z1/2  R 
1  2  12 


(D-31 


If  it  is  assumed  that  the  correlation  coefficient  Is  approximately 
unity  in  the  interaction  region  of  the  parametric  receiver,  and  if 
expressions  similar  to  Eq.  (D-31)  are  written  for  the  terms  <s+1sp2>  and 
<Sx1Sio>»  then  Eq.  (D-30)  becomes 


■4 


oo 


<SpE>  i  <u2>  k2  I  R(p)dp 


M/[ 


zj/2  z^2  +  zJ/2(L-z2)1/2  +  (L-z/'2  Z^2 


+  (L-Zl)1/2  (L-z„) 


1/2 


dz^dz2 


Integration  is  straightforward  and  yields 


OO 

<S2r>  =  1.778  <u2>  klLJ R(p)dp 


Approximating  the  correlation  term  R(p)  by  the  Gaussian  function 
2  2 

exp(-p  /a  )  gives 


=  0.8889  /n  <y 


2 


> 


k2  aL 
P 


(D-32) 


This  result  relates  the  mean  squared  phase  fluctuations  in  the  output  of 

2 

a  mobile  parametric  receiver  to  the  medium  parameters  <p  >  and  a. 
Equation  (D-32)  is  valid  only  for  large  values  of  the  wave  parameter 
(D>>1),  as  indicated  in  Eq.  (C-l) .  For  small  values  of  the  wave  param¬ 
eter  (D<<1),  the  phase  fluctuations  given  by  Eq.  (C-l)  increase  by  a 

35  2 

factor  of  two,  and  the  result  for  <SDD>  becomes 

rK 


> 


1.778  Sn  <\i2>  k2  aL 
P 


(D-33) 
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G. 


Summary  and  Discussion 

Theoretical  analysis  of  parametric  reception  with  turbulence  present 
between  pump  and  hydrophone  has  produced  the  following  expressions  for 
mean  squared  amplitude  and  phase  fluctuations. 

Amplitude  Fluctuations 


> 


„  1MC  p2  ,7  6  T 11/6 
0.1415  C  k  '  L 
n  p 


Phase  Fluctuations  for  D>>1 


> 


0.8889  /if  <y2>  k2  aL 
P 


(D-29) 


(D-32) 


Phase  Fluctuations  for  D<<1 

<S^  >  =  1.778  /if  <y2>  k2  aL  (D-33) 

r  K  P 

A  number  of  assumptions,  made  in  deriving  these  results,  are 
summarized  in  Table  D-l.  Two  assumptions,  (6)  and  (8),  are  particularly 
important  in  determining  the  validity  of  the  theoretical  results.  Assump¬ 
tion  (6),  which  states  that  the  fluctuations  are  small,  effectively  limits 
the  levels  of  rms  amplitude  or  phase  fluctuations  that  can  be  predicted 
to  approximately  10%  of  the  mean  amplitude  or  phase,  respectively;  i.e., 


> 


1/2 


<  0.1 


and 


1/2 


<  0.1 
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TABLE  D-l 

LIST  OF  ASSUMPTIONS 


1.  Spherical  pump  wave  and  planar  signal  wave. 

2.  Frequency  broadening  of  waves  due  to  scattering  is  sufficiently 

small  that  a  ,  a  ,  and  a.  are  constants, 
s  p  ± 

3.  Amplitude  and  phase  fluctuations  of  the  signal  wave  are  negligibly 
small. 

4.  There  is  complete  correlation  Bp,  Sp,  B+,  and  S+  along  the  spherical 
wave  front  in  the  interaction  region. 

5.  Sphericity  of  pump  wave  is  imall  enough  to  assume  "cophasal  source 
wafers"  in  the  interaction  region. 

6.  The  scattering  is  sufficiently  weak  that  BpR,  SpR  B  ,  S  ,  B+,  and  S+ 

remain  small  compared  to  unity.  >  P  P  - 

7.  a  +a  -a  =0. 

s  p  ± 

8.  Fluctuations  Bp,  Sp,  B+,  and  S+  have  complete  longitudinal 

correlation  in  the  interaction" region  (R  =R  ,=R  =1). 

P  P,±  r 

Turbulence  in  the  interaction  region  is  homogeneous  and  isotropic. 


9. 


In  some  cases  of  strong  turbulence,  high  pump  frequency,  or  long  array 
lengths,  assumption  (6)  may  be  violated,  so  that  predicted  results  become 
inaccurate.  A  simple  method  of  extending  the  range  of  validity  for  the 
theory  is  discussed  in  Appendix  E. 

Assumption  (8),  which  states  that  the  pump  and  interaction  frequency 

fluctuations  have  complete  longitudinal  correlation  in  the  interaction 

region,  has  great  impact  on  the  accuracy  of  the  theoretical  results.  It 

can  be  seen  from  Eqs.  (D-24)  and  (D-31)  that,  if  these  correlation  terms 
2  2 

are  zero,  then  <B_n>  and  <St)_>  become  zero.  If  the  correlation  terms  are 
PR  PR 

unity,  as  assumed,  then  the  results  given  above  [Eqs.  (D-29) ,  (D-32) ,  and 

(D— 33) ]  are  valid.  Thus  the  results  obtained  in  the  previous  two  sections 

are  worst-case  results.  If  the  pump  and  interaction  frequency  fluctuations 

are  to  some  degree  uncorrelated  over  the  interaction  region,  the  fluctua- 
2  2 

tions  <B_._>  and  <S__>  are  reduced,  and  must  be  calculated  using  the 
PR  PK 

general  expressions  in  Eqs.  (D-24)  and  (D-30) .  The  correlation  coeffi¬ 
cients  R  ,  R  +,  and  Rr  are  therefore  very  important  in  determining  the 

value  of  <B  >  and  <SDti>.  Values  of  these  correlation  coefficients  depend 
PR  PR 

strongly  upon  the  intensity  and  geometry  of  the  turbulence  present  in  the 
interaction  region,  and  will  require  experimental  results  for  accurate 
determinations . 
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APPENDIX  E 

EXTENSION  OF  RESULTS  FOR  STRONGER  TURBULENCE 


2  2 

A  method  for  obtaining  approximate  values  of  <B  >  and  <Sno>  for 

*  R  PR 

fluctuations  greater  than  10%  of  mean  value  is  discussed  in  this  appendix. 

Comparison  of  Eqs.  (D-27)  and  (D-29)  shows  that,  for  the  case  of 

weak  fluctuations  (B  <<1  and  B  <<1),  the  mean  squared  amplitude  fluctua- 

p  i  K 

tions  in  the  pump  and  interaction  frequency  waves  are  approximately  equal. 
In  other  words. 


> 


<BPR> 


(E-l) 


according  to  the  results  of  the  analysis  in  Appendix  D.  Similarly, 
using  Eqs.  (C-2)  and  (D-32) , 


<s2>  =  <s2  > 

p  PR 


(E-2) 


although  in  this  case  the  approximation  is  much  cruder  than  that  in 
Eq.  (E-l). 

Equations  (E-l)  and  (E-2)  form  the  basis  of  a  simple  theoretical 
model  that  can  be  used  to  make  approximate  predictions  of  fluctuations 
when  the  small  perturbation  results  presented  in  Appendix  D  are  no  longer 
valid.  It  is  assumed  that  the  approximations  in  Eqs.  (E-l)  and  (E-2) 
remain  valid  regardless  of  the  intensity  of  turbulence  in  the  interaction 
region.  It  should  be  noted,  however,  that  Eqs.  (D-29)  and  (D-32)  assume 
that  the  pump  fluctuations  are  correlated  throughout  the  interaction 
region.  If  this  assumption  is  violated,  then  approximations  (E-l)  and 
(E-2)  become  less  accurate. 

2 

When  Eqs.  (E-l)  and  (E-2)  are  valid,  then  the  determination  of  <B  > 

2  FK 
and  <S  >  reduces  to  a  problem  of  determining  the  amplitude  and  phase  flue- 
rR 

tuations  of  a  linear  acoustic  wave  (the  pump  wave)  propagating  in  a  tur¬ 
bulent  medium.  A  solution  to  this  problem  has  been  obtained  by  the 
method  of  smooth  perturbation,  or  Rytov's  method,  and  is  available  in  the 
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t 


26  27  31 

literature.  ’  *  Results  for  the  amplitude  and  phase  fluctuations 

are 


> 


<B2>  =  0.13  C2  k776  L11^6 
P  n  p 


(E-3) 


and 


<S2  >  =  <S2>  =  0.50  /?  <u2>  k2  aL.  (E-4) 

PR  p  p  ' 

Major  assumptions  used  in  deriving  Eq.  (E-3)  are  that  the  turbulence  is 
isotropic  and  homogeneous,  that  the  pump  wave  is  spherically  spreading, 
and  that  only  the  inertial  subrange  of  the  turbulence  spectrum  contributes 
to  the  acoustic  fluctuations  [see  the  discussion  related  to  Eqs.  (C— 7 ) — 
(C-10)  in  Appendix  C].  Equation  (E-4)  was  derived  assuming  that  the  cor¬ 
relation  coefficient  of  the  refractive  index  variations  is  Gaussian  and 
that  the  wave  parameter  D  is  much  greater  than  unity  [see  Eq.  (C-2)]. 

Rytov's  method  is  valid  for  rms  fluctuation  levels  up  to  about  0.5; 
i.e.,  the  fluctuations  must  be  such  that 


<BPR> 


1/2 


<  0.5 


and 


> 


1/2 


<  0.5 
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TURBULENCE  PARAMETERS  FOR  THEORETICAL  EXAMPLES 
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In  this  appendix  expressions  are  derived  for  the  turbulence 

2  1/2  2  1/2 

parameters  that  are  used  in  calculating  values  of  <B_„>  and  <SD_>  . 

rK  PR 

The  structure  constant  C  in  terms  of  the  refractive  index  variations 

n  43 

and  the  outer  scale  of  turbulence  is  given  by 

2  2  -2/3 

-  1.91  <v  >  Lo^J  ,  (F-l) 


where 


Ac 


Ac  is  the  local  variation  in  sound  speed,  and 

L  is  the  outer  scale  of  turbulence,  determined  by  the  dimensions 
of  the  flow. 

In  estimating  a  value  of  y  for  the  turbulent  wake  of  the  pump  transducer, 

44 

the  results  given  by  Hinze  for  the  wake  of  a  cylinder  may  be  used  as  an 
approximation.  For  the  central  part  of  the  wake,  the  normalized  mean 
squared  velocity  fluctuations  in  the  direction  of  mean  flow  are 


<(Av)2> 


=  0.1 


(F-2) 


where  v  is  the  maximum  velocity  behind  the  obstacle  in  the  direction  of 
flow.  Rearranging  Eq.  (F-2)  gives 


<(Av)2> 


1/2 


=  /U  v  =  0.3162  v 


(F— 3) 


Because  Av  is  in  the  direction  of  the  pump  wave  propagation,  the  local 
sound  speed  is  increased  by  an  amount 


<(Ac)2>  *  <(Av)2> 


(F-4) 


assuming  that  convection  is  the  only  mechanism  contributing  to  a  change 
of  sound  speed  in  the  wake.  Using  Eqs.  (F-3)  and  (F-4),  the  mean  squared 
refractive  index  variations  are 


(F-5) 


<jj2>  - 


<(Av)^> 


If  Che  mean  sound  speed  is  approximately  equal  to  1500  m/sec ,  then  using 
Eqs.  (F-l)  and  (F-5),  the  structure  constant  may  be  written  as 

C2  =  8.487  x  10-8  v2  L~2/3 
n  o 

Another  variable  that  needs  to  be  expressed  in  terms  of  the  wake 
parameters  is  the  refractive  index  correlation  length  a.  One  simple  way 
is  to  assume  that 


k  =  0.5(k  k  ) 
t  pm 


1/2 


2tt 

a 


(F-6) 


This  may  be  interpreted  physically  as  setting  the  correlation  length 
equal  to  the  eddy  sizes  that  have  maximum  effect  upon  the  sound  propaga¬ 
tion.  This  approximation  is  rather  crude,  but  may  be  tested  when  meas¬ 
urements  of  turbulent  flows  are  made.  If  the  wave  number  ic  is  related 

m 

to  the  outer  scale  of  turbulence  by 


2tt  4tt 


2 

then  a  may  be  found  from  Eqs.  (F-6)  and  (F-7)  to  be 

a  =  4.737  x  io-2  ST 
o 


(F-7) 


where 


40 


k  =  5600  m-1  at  30°C 
P 


t 
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